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ABSTRACT 

Analysis of the peculiar velocities of galaxies should take account of the uncertainties in both redshifts 
and distances. We show how this can be done by a numerical application of the action principle. The 
method is applied to an improved catalog of the galaxies and tight systems of galaxies within ih^ Mpc, 
supplemented with a coarser sample of the major concentrations at Mpc to 20h^ Mpc distance. 

Inclusion of this outer zone improves the fit of the mass tracers in the inner zone to their measured 
redshifts and distances, yielding best fits with reduced x 2 m redshift and distance in the range 1.5 to 2. 
These solutions are based on the assumption that the galaxies in and near the Local Group trace the 
mass, and a powerful test would be provided by observations of proper motions of the nearby galaxies. 
Predicted transverse galactocentric velocities of some of the nearby galaxies are confined to rather narrow 
ranges of values, and are on the order of 100 km s _1 , large enough to be detected and tested by the 
proposed SIM and GAIA satellite missions. 

Subject headings: cosmology: theory — galaxies: distances and redshifts — Local Group 



1. INTRODUCTION 

This paper continues the analysis of the dynamics of rel- 
ative motions of the galaxies in and near the Local Group. 
In previous numerical applications of the action principle 
either the present distances of the assumed mass tracers 
are given, and the measure of goodness of fit is the differ- 
ence between model and catalog redshifts ( Peebles- 1989), 
or the present redshifts are given, and the measure of good- 
ness of fit is the difference between model and catalog dis- 
tances]^] Since errors in redshift and distance both are 
important we have adapted the numerical action method 
to allow relaxation to a stationary point of the action and 
a minimum of the sum of the squares of the weighted dif- 
ferences of model and catalog redshifts and distances. The 
method is described in § 3. 

In § 4 we present an application of the method to the 
catalog discussed in § 2. Recent advances in measure- 
ments of distances of nearby galaxies are incorporated in 
a catalog of the distances and radial velocities of the galax- 
ies and tight concentrations of galaxies within ihj^ Mpc. 
(Hubble's constant is H Q = 75/175 km s _1 Mpc ). This 
inner zone catalog is supplemented by a coarser catalog 
of distances and radial velocities of the main galaxy con- 
centrations in an outer zone at 4~20h^ Mpc distance. 
We work with the assumption that the luminosities of the 
galaxies trace the mass. Two control cases support this as- 
sumption: (i) inclusion of the outer zone in the dynamical 
solutions improves the fit to the measured redshifts and 
distances in the inner zone, and (ii) replacing the angular 
positions by random positions on the sky, while keeping 
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the catalog redshifts and distances, worsens the fit. 

The application presented here is meant to illustrate the 
new numerical method and catalogs; we hope to present a 
systematic study of parameters in a later paper. With the 
parameter choices for cosmology and mass-to-light ratios 
(M/L, using blue-band light) listed in § 4 we get values of 
the reduced x 2 for distances and redshifts that in the best 
cases are between 1.5 and 2, close enough to unity to add 
to the evidence that our model is a useful approximation to 
reality. A considerably stronger test, from measurements 
of transverse motions of nearby galaxies, is discussed in 
§5. 

2. CATALOGS OF NEARBY GALAXIES 

Our compilation of mass tracers and potential targets 
for measurements of proper motions is restricted by three 
considerations. (1) Accurate distances are needed to con- 
strain models. (2) In this first exploration of the method 
we want to restrict the number of orbits to reconstruct. (3) 
Complex — strongly nonlinear — orbits can only be re- 
covered with an extensive exploration of initial conditions. 
The first two led us to study about 40 objects. The third 
compels the merging of some closely adjacent galaxies into 
single entries in our catalog. 

Regarding the last point, it is possible in principle to find 
numerical action solutions in cases where galaxies have al- 
ready made close passages, as long as there has not been 
significant exchange of orbital and internal energy. For 
example, the interactions between M81, M82, and NGC 
3077 left extended tidal streams of neutral hydrogen that 



provide d etailed constraints on the relative motions (Yun. 
~Tc] 1994), a case that would be fascinating to study. 
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However, the current work has a much less ambitious goal. 
If galaxies are so near to each other that the orbits may 
be complex, we merge the information on these systems: 
the catalog entries are the sums of the luminosities and 
the luminosity- weighted positions and velocities. What is 
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meanl by 'near to each other' depends on the available in- 
formation. For galaxies that are close to us we have better 
relative discrimination of positions and a more complete 
census of the major mass influences. Consequently, we at- 
tempt to explore somewhat more complex situations in the 
nearest volume. 

Our catalog has 'inner' and 'outer' zones. The former, 
within Mpc distance, is the main sample in this 

study. About 90% of the light in the inner zone resides 
in just 8 of the objects. If mass and light are strongly 
correlated then most of the inner catalog mass is asso- 
ciated with these few giant systems. It is important for 
the dynamical analysis that these few big objects be lo- 
cated as precisely as possible in position and velocity. The 
many small galaxies make a minor contribution to the light 
and, we are assuming, to the mass. We seek to include as 
many of these objects as feasible because they are valuable 
probes of the gravitational potential, but their exclusion 
would have little effect on the dynamics. Thus we include 
all the massive galaxies in the inner zone, whatever the re- 
liability of their distance estimates, but only dwarfs that 
are sufficiently isolated to have simple orbits and for which 
we have good distance estimates. 

Thanks especially to the Hubble Space Telescope, rea- 
sonably good distances are emerging for many galaxies in 
and at the fringe of the Local Group. Two techniques 
are providing most of the data: the correlation between 
pulsa tion period and luminosity of Cepheid variable stars 



fe Fisher (1987) provide maps of these structures. 



(eg., Freedman ct al. 1994) and the calibrated luminos- 
ity of the tip of the Red Giant Branch for evolved low- 
metallicity stars (eg., Cole et al. 1999). Each can provide 
distances with 10% standard deviation. There is a fur- 
ther 10% uncertainty in the underlying zero points, but 
this normalization is close to a common factor in all dis- 
tances and so enters only as a scale factor in the dynam- 
ics. The distance scale used in this paper is con sistent 
with H = 75 km s" 1 Mpc" 1 ( [Tully fc Picrec] 2000; |Tonry 
et al{] 2000); a zero point shift is accommodated by the 
factor /175 = Ho/75. Distance measurement coverage of 
individual galaxies falls rapidly beyond Mpc. 

Some comments on the characteristics of the nearby 
structure are in order. Almost all galaxies within lhj^ Mpc 
have negative systemic velocities, which suggests this re- 
gion is gravitationally bound but not yet relaxed. The 
galaxies in this region make up the Local Group. This con- 
densation is part of an expanding filamentary structure, 
with the nearest big galaxies close to each other in the sky 
in what are called the Maffei-IC342 and M81 groups. The 
filament originates in the northern Galactic hemisphere at 
a conjunction of filaments near the Ursa Major and Coma 
I clusters, passes through small knots of galaxies in the 
constellations of Canes Venatici and Ursa Major, through 
our position, and can be followed in the southern Galactic 
hemisphere out to the NGC 1023 Group. The thread of 
galaxies that makes up what has been called the Sculptor 
Group can be viewed as a minor fragment of the main fil- 
ament. It seems this filament is paralleled by a second one 
that originates in the Galactic north at the Virgo Cluster, 
passes through a big knot including the Sombrero galaxy, 
then passes nearest to us at the Centaurus Group, through 
the Galactic plane in the vicinity of Circinus galax y, and 
onward to meet with the Telescopium-Grus Cloud. Lully 



The outline of these features can be seen in Figure 1, 
though in skeleton form because we have plotted only the 
important dynamical constituents and those that serve as 
useful test particles. Our inner zone catalog in Table 1 
is shown as the filled symbols. This sample is thought to 
be complete in luminous objects (Lb > 10 10 Lq) within 
4/iyg Mpc, and we are assuming it is a good sample of 
the mass. The unnamed objects are dwarfs with decent 
distance measurements. All the inner zone objects are 
concentrated near the equatorial plane in supergalactic 
coordinates. If light traces mass the major gravitational 
influence on the Local Group from the inner zone is the 
Maffei/IC342/M81 complex, with some influence from the 
Centaurus Group, a minor effect from the modest galax- 
ies in the Sculptor Group, and little else. There are no 
important mass elements near the boundary of the inner 
zone. This degree of isolation and the sharp decline in the 
completeness of distance measurements led us to place the 
boundary for the main inner zone at 4ft," 1 Mpc. 

Gravitational interactions with more distant mass con- 
centrations affect the dynamics within our inner zone. The 
nearest substantial exterior concentrations of galaxies are 
at distances of 15-20/ifg 1 Mpc. Then at 30-60/ifg 1 Mpc 
there arc still larger concentrations in the Norma-Hydra- 
Centaurus and Perseus-Pisces supercluster regions. For 
this paper we ignore the big structures beyond 30h^ Mpc 
but try to account for the significant structures within an 
outer zone at A-20h^ Mpc distance. The distribution 
of important galaxy concentrations in this outer zone is 
simple enough to be modeled by the 14 objects in Ta- 
ble 2. These groups or clusters (numerical names from 



Tullyj 1988) are shown as open symbols in Figure 1, and 



the eight that cause the largest gravitational tides at our 
position (based on mass divided by distance cubed) are 
identified by name in the bottom pair of panels in the fig- 
ure. 

The Virgo Cluster is not an outstandingly dominant 
contributor to the local luminosity, but there is good ev- 
idence that it is a larger contributor to the mass. Nu- 
merical action modeling, which mostly refers to the low 
density environments where most galaxies are found, in- 
dicates the cosmological mass density parameter is fl m = 
0.25 ± 0.2 (95% formal probable error) with mass-to-light 
ratio M/L ~ 150^7^ M^/Lp, in typic al gal axy environ- 
ments ( |Shaya, Peebles, fc Tull^ 1995; |Tullyj 1999). The 



evidence that there is much more mass in the Virgo Clus- 
ter than would be inferred from this value of Mj L includes 
the velocity dispersions of the early and late-type galaxies 
near the core, the high veloc ities of galaxies infalling on 
first approach to the cluster ( Tully fc Shaya 1984; 1998), 
and the virgocentric flow derived from the distances in the 
surface brightness fluctuations survey ( Tonry et al. 2000). 



Tully fc Shaya (1998) conjecture that the Virgo Cluster 



is not the only nearby place where M/L is large, that re- 
gions with short dynamical times (t co u apse ~ 10 9 years) 
have high M/L, and that such regions are readily identi- 
fied as nests of elliptical and SO galaxies. It is argued that 
these regions have M/L at least five times that of regions 
with long dynamical times {t co u apse ~ 10 10 years). If so, 
four other systems with large M/L within 20h^ Mpc are 
the E/S0 knots associated with the Fornax, Dorado, Coma 
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Fig. 1. — Projected locations of catalog objects. Top pair of 
panels: Region within ±4/i^" 8 Mpc from two orthogonal viewpoints 
in supergalactic coordinates. Bottom pair of panels: Region within 
±20/if 5 1 Mpc from the same orthogonal viewpoints. Shaded wedges 
identify zones of Galactic obscuration. 



I, and Leo clusters/groups. These entities are located by 
large boxes in the bottom two panels of Figure 1; we as- 
sign them M/L five times that of the nine low density 
groups with long dynamical times and mainly late-type 
galaxies in the outer zone. The most massive of the low 
M/L members of the outer zone, the Sombrero and Ursa 
Major clusters, are identified by the larger open circles in 
Figure 1. The NGC 6946 group is the nearest of the outer 
zone systems at 6h^ Mpc (a distance that still is quite 
uncertain). 

The mass we assign to the Virgo Cluster is four times 



that of the second most massive system in Table 2, and 
larger than the sum of all the other masses within 20hj^ Mpc. 
A simple measure of the interactions among inner and 
outer zone objects is given in the last column in Table 2: 
the relative contribution of each outer object to the present 
tidal field at our position. One sees again the importance 
of the Virgo cluster, and the not insignificant sum of con- 
tributions from the rest of the outer zone objects. 

The luminosities of all the galaxies in the outer zone 
that are not in the 14 objects in Table 2 add up to a 
luminosi ty density that is compar able to the cosmic mean 
used by |Shaya, Peebles, fc Tullyj (1995). That is, the 14 
mass tracers in the outer zone approximate the overdensity 
of luminosity, and plausibly of mass, within 20h^ Mpc. 

Our numerical analysis assumes distance measurement 
errors are symmetrically distributed in the distance mod- 
uli, with the standard deviations o~ u listed in Tables 1 
and 2. The standard deviations a cz in the redshifts in 
Table 1 are measurement errors, while the a cz listed in Ta- 
ble 2 are mean deviations of the velocity dispersions within 
these clusters and groups of galaxies. A crude model for 
the motion of each galaxy relative to the dark halo it is 
supposed to trace is presented in § 4. 

The objects near the outer boundary of the outer zone 
likely are strongly influenced by structures beyond 20h^ Mpc 
that are not in our analysis. Hence although orbit com- 
putations arc constrained by the distances and velocities, 
with their uncertainties, in both zones, the solutions are 
evaluated only from the fits to the redshifts and distances 
in the inner zone. 

3. SOLUTIONS THAT MINIMIZE x 2 

3.1. Numerical Action Method 

This initial discussion, for simplicity, deals with a single 
particle. A solution to the equation of motion generally 
has neighboring solutions with neighboring values of the 
present distance and redshift, and we seek the member of 
this family that minimizes 



X 2 = (z- z c ) 2 /a 2 + - He) 2 1^1 



(1) 



where the redshift and distance modulus in the solution are 
z and h, the measured catalog values are z c and /i c , and 
their standard deviations are o~ z and <7 M . The relaxation of 
a single-particle orbit to stationary points of the action and 
of x 2 is applied separately to each member of the catalog 
and iterated until the gradients of the action and x 2 vanish 
to machine accuracy. Relaxation to a minimum of the sum 
of the x 2 over all particles (except the reference), rather 
than iterated relaxation to minima of single-particle x 2 , is 
feasible but the matrix inversion is a considerably heavier 
computation. 

We represent the orbit of a particle by its comoving 
positions x„ at a sequence of time steps t n , 1 < n < n x , 
together with the present distance d at t ime t Q — t nx +x 
and given present angular position, as in [Peebles] (1995). 
The action function is 



A = S + H a d 2 /2 — d(cz + vmw • x), 



(2) 



where 



71=1 



a n+l/2 (Xn+i - X n ) 2 
{t n +l ~ tn) 
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{t — ^+1/2)^0 



(3) 



The term S approximates the time integral of the kinetic 
minus the potential energies. The potential for particle i 
is 

v = -e E^-^- W 



4a 



The physical distance between particles i and j is a(t)xij. 
The second term in refers the potential to the back- 
ground cosmological model with present Hubble constant 
H and density parameter Q. in matter with mass density 
that varies inversely as the cube of the expansion param- 
eter a(t). The terms added to S in equation (||) represent 
a canonical transformation that in effect exchanges radial 
positions and momenta, thereby changing the boundary 
condition from fixed present distance to fixed pres ent ra- 



dial velocity (the details of which arc presented in Phclpi 
2000). The present radial velocity of the particle relative 
to the Milky Way is cz, and the last term in equation (J2J) 
adds the component of the peculiar velocity vmw of the 
Milky Way in the direction x of the particle. (The veloc- 
ities are such tha t the center of mass moves according to 
eq. [7] in Peebles 1994). The derivative of the action with 



respect to the present distance is 



dA _ a nx +i/2 x ■ (x - x„, 
dd a a t Q - t nx 

(t t \ dV ° 

+ H Q d — cz — vmw • x. 



(5) 



The first term on the right hand side is the radial canon- 
ical momentum of the particle a half time step before the 
present, where x is the unit vector to the present particle 
position. The second term brings the radial momentum 
to the present epoch. When dA/dd = this agrees with 
the radial peculiar velocity given by the last three terms. 
One similarly sees that the conditions dA/dx n = rep- 
resent th e equation of motion in leapfrog approximation 
( |Peeblesj 1995). The motion of the Milky Way is at a sta- 
tionary point of the function S of the coordinates of this 
orbit, with the boundary condition that the Milky Way 
ends up at the origin. 

To simplify discussion of the relaxation to stationary 
points of A and x 2 let x a represent the 3n x + 1 coordinates 
for the orbit of a particle (or the 371^ coordinates for the 
Milky Way), with a = for the present distance xq = d. 
The first and second derivatives of the action with respect 
to these variables are A a and A a p. The square of the 
gradient of the action for particle i is I\ = ^2(A a ) 2 . The 
coordinate shift 



8x a 
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(6) 



reduces Tj if e is small enough. If x a is close to a stationary 
point, so A is close to a quadratic function of the distance 
from the stationary point, then e = 1 brings the orbit much 



closer to the stationary point, whether a local extremum 
or saddle point . 

Now consider neighboring stationary points of the ac- 
tion, where A a — 0, belonging to neighboring redshifts z 
and z + dz. The result of differentiating A a = with re- 
spect to z, and remembering that the coordinates 
functions of redshift and that z also appears explicitly in 
A (eq. |) is 

AolP^T- = cS a .Q. (7) 

dz 

We see that in a family of solutions related by continuously 
varying the redshift z the derivative of the present distance 
d with respect to the redshift is 



dd 

dz 



cA~ 



00- 



(8) 



The inverse of the matrix of second derivatives of A thus 
shows how to move the orbit toward a stationary point of 
A (eq. ||) and how the present distance at a stationary 
point of A changes when the redshift is adjusted. We use 
the latter to relax to a minimum of % 2 . The derivative of 
X 2 (eq. 0) with respect to the redshift satisfies 



where 



^ 2 9x2 10 
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dd In d/d c 
dz d 



(9) 



(10) 



We find that it is a good numerical approximation to sup- 
pose dd/dz is independent of z. In this case the second 
derivative satisfies 



<9 2 y 2 



dz 



R 



dd 

dz 



(l-lnd/4). (11) 



The analog of equation (pi) for a redshift adjustment that 
moves the orbit toward the minimum of x 2 is 



Sz 



dx 2 /dz 
d 2 x 2 /dz 2 



(12) 



To sample the families of solutions we start from ran- 
dom orbits with the positions x„ at 1 < n < n x placed 
independently at random in the sphere that contains the 
mass of the system of objects in the background cosmolog- 
ical model. The iterative relaxation first uses the matrix 
inverse A~ 1 a p to find dd/dz and hence the redshift ad- 
justment Sz (and the accompanying adjustment to dA/dd) 
that moves the orbit toward a minimum of x 2 , and then 
uses the same matrix inverse to find the coordinate adjust- 
ments 5x a that move the orbit toward a stationary point 
of A, 

3.2. Numerical Fixes 

We expect colleagues who wish to explore applications 
of this method would do well to seek independent ways 
to address the numerical problems arising, so as to im- 
prove our inelegant and likely inefficient fixes, but to aid 
reproducing our results we list the main elements of our 
procedure. 
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The square of the gradient of the action summed over 
all objects i and position coordinates a is 



r = 53r i = 53(5A/fl a 



(13) 



With the velocity unit 100 km s _1 , length unit 1 Mpc, 
and mass unit 1 x 10 11 solar masses, the square of the 
gradient of the action for the initial random orbits typi- 
cally is T ~ 10 8 . At T < 10 -10 further iterations usu- 
ally reduce the square of the gradient of A to zero to ma- 
chine accuracy (32 bit) without appreciable perturbations 
to the orbits. In the solutions presented here the iteration 



stops at T ~ 10~ zu , at which point (dx i ' /dz) z is similarly 
small . The iteration can approach a point where l l lias 
a local minimum but not a zero. In our experience when 
this happens T is dominated by the contribution from one 
object, so when Ti > T/2 for object i its orbit alone is 
adjusted until Tj < 0.01T or to a maximum of 300 itera- 
tions. If this happens for particle i more than 100 times 
we choose a new random orbit for the object and relax 
it alone until Ti < 0.01T. Negative distances, and orbits 
that cross at zero separation, are unphysical but allowed 
by the mathematics. During the relaxation a distance can 
move from negative to positive, so when object i has dis- 
tance di < 30 kpc and r, < 1 we choose a new random 
orbit for this object and relax it until Tj < 0.01T. When 
r < 1 and the minimum comoving separation of a pair of 
particles is less than 30 kpc we apply the above procedure 
to the less massive one. Equation (||) assumes the coor- 
dinates are at a stationary point of A, but we apply the 
redshift adjustment at each coordinate adjustment start- 
ing from random orbits. When the coordinates are far 
from a stationary point the approximation to d 2 \ 2 /dz 2 in 
equation ( |l"ij ) can be negative. In this case we relax down 
the gradient of % 2 , using 



Sz = —( 



d\ 2 /dz 1 



(14) 



When the second derivative of \ 2 is positive we use the 
smaller of the redshift adjustments from equations (12) 
and (H|). There are many solutions to the equation of 
motion for the objects in Tables 1 and 2; we seek those 
with rcdshifts and distances that are close to the catalog 
values. When the relaxation brings T below 10~ 5 each ob- 
ject in turn with individual \ 2 greater than 25 is given 40 
attempts at a new orbit, relaxed from a random one, to 
bring the individual \ 2 below 25. After this treatment of 
each offending orbit all orbits are iteratively relaxed un- 
til r < 10~ 5 again. The operation is repeated, allowing 
a maximum of 5 sets of 40 attempts for each object. In 
the analysis of the real catalogs usually all individual x 2 
are below 25 after about 1000 iterations through all or- 
bits. In the control case with randomly placed present 
angular positions the allowed number of attempts usually 
is exhausted by 1000 iterations, leaving some objects with 
quite large \ 2 - At n x > 10 time steps the computation 
time scales as n%, dominated b y the matrix inversion codes 
LUDCMP and LUBKSB from press et al.| (1992). 



4. APPLICATION TO THE NEARBY GALAXIES 
4.1. Parameters and Control Samples 



In this illustration of the numerical method we fix 
some parameters by considerations other than the mo- 
tions of the nearby galaxies. We adopt a cosmologically 
flat Friedmann-Lemaitre model with matter density and 
Hubble parameters 



H m = 0.25, 



Ho = 75/i75 km s -1 Mpc 1 . 



(15) 



Our choice of £l m is near the lower end of the range de- 
rived from recent measurements of the angular fluctuations 



in the thermal cosmic background radiation (eg. Hu et al 
2000) and near the upper end from analy ses of gal axy dy 



namic s (as summarized by Tully & Shaya 1998 and Bahcalj 



et al 
of ga 



2000). The value of H fits the recession velocities 
axies at cz = 5000 to 8000 km s _1 that plausibly are 
in the Hubble flow ( [Tully fc Pierce] 2000). The scale factor 
/175 takes account of a close to common uncertainty in the 
distance zero points for H a and the distances in Tables 1 
and 2. 

Luminosities are computed at the catalog distances. 
The masses in the inner zone in Table 1 assume the mass- 
to-light ratio 



M/L(inner) = 75h 75 M G /L Q . 



(16) 



This choice makes the sum of the masses of M31 and 
the Milky Way consistent with the spherical model for 
their relative motion, and the total mass contrast at 
r < Ahjs Mpc is SM/(M) = 1.0, consistent with the mod- 
est gravitational slowing of the local expansion rate just 
outside the Local Group. Following the discussion in § 2, 
we take account of the light in the outer zone that is not 
in the mass concentrations in Table 2 by adopting a larger 
mass-to-light ratio, 



M/L(outer) = W0h 75 M G /L 



(17) 



for the late-type systems in Table 2, and 75O/175 Mq/Lq 
for the early type systems. This makes the contrast 
SM/(M) = 0.6 at r < 20h^ Mpc, again a reasonable 
number. 

To take account of the distributed mass around the trac- 
ers we change the first term in the gravitational potential 
in equation ([!]) to V oc (xi 2 + (cj + a,) 2 ) -1 / 2 , with comov- 
ing cutoff length Cj = 30 kpc for all particles in the inner 
zone, 2 Mpc for Virgo, and 700 kpc for all other objects 
in the outer zone. The cutoff is larger in the outer zone 
because these objects are meant to represent more broadly 
distributed mass. When all cutoff lengths are set to zero it 
allows more sharply curved and perhaps unphysical orbits, 
but the statistics are little changed. 

The standard deviations er cz in redshift take account of 
the dispersions of velocities in the large systems of galaxies 
in Table 2, but in Table 1 reflect only the measurement 
errors. Since we can only guess at the typical difference 
between the velocity of the galaxy and the motion of the 
center of mass of the dark halo it is supposed to represent, 
we write the standard deviation in redshift as 



= K, + O 



2U/2 



(18) 



and we compare results for a — and ca Q = 15 km s _1 . 

We label dynamical solutions for the combined set of 
orbits of the objects in both zones as cO and cl5, for the 
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two values of a . In the first set of control cases, labeled iO 
and il5, the outer zone has been removed from the dynam- 
ics. In the second set of control cases, labeled rO and rl5, 
the outer zone is included and the catalog distances and 
redshifts are used but each object is placed at a randomly 
chosen position in the sky (with the same angular position 
in all solutions). All statistics for all cases are based on 
the objects in the inner zone alone. 

We use n x = 20 time steps uniformly spaced in the ex- 
pansion parameter a(t) and e = 0.2 in the relaxation equa- 
tions @, (|l2|), and (|l4|). The statistics are quite similar 
at e = 0.1 and e = 0.3, and at n x = 15 and n x = 30. We 
compute 30 solutions from random orbits for each case. 
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reduced x z 



Fig. 



2. — Distributions of reduced \ 2 m distance and redshift 
for 30 solutions for each of the six cases. In the broken histograms 
cr = 0, in the solid, ca a = 15 km s _1 . 



4.2. Numerical Solutions 

Figure 2 shows the distributions of reduced % 2 in redshift 
and distance summed over the 29 objects in the inner zone 
(excluding the Milky Way) for the 30 solutions for each 
of the six cases. Increasing the standard deviations in the 
redshifts makes it considerably easier to find solutions with 
reduced x 2 l ess than 3. Removing the outer zone from 
the dynamics or scrambling the angular positions makes it 
more difficult to find these relatively good fits to the data. 
The differences of the best \ 2 among the 30 solutions for 
each case are much smaller than the differences among 
the distributions of x 2 from all 30 solutions, but the best 
solution for the real catalogs is consistently better than 
the best solution for the control cases. 

The last two columns in Table 1 and the filled squares in 
Figure 3 show the normalized differences between model 
and catalog redshifts and distances for the 29 objects in 
the best solution for the real catalogs with a Q = (case 
cO). For comparison we show as the open symbols the nor- 
malized differences in the best of the solutions for scram- 
bled angular positions (case rO). The normalized scatter is 
smaller in redshift and correlated with the scatter in dis- 




FlG. 3. — Distributions of normalized differences between model 
and catalog redshifts and distance moduli in the best solutions for 
cases cO (filled squares) and rO (open squares). 



tance. This is because dd/dz (eq. ||), the rate of change 
of redshift with respect to distance in a family of solutions, 
usually is positive and, in the units in Figure 3, the slope 
usually is greater than unity. The single-particle % 2 are 
minimized at the closest approach of the trajectory of red- 
shift as a function of distance to the origin in Figure 3, 
so the minima scatter around a line that trends down to 
the right. Some trajectories happen to pass close to the 
origin in Figure 3; their relaxation to minimize x 2 pro- 
duces unrealistically good fits to the catalog redshifts and 
distances. 

Figure 4 shows the distributions of single-particle % 2 
in redshift and distance for the best solutions for each of 
the six cases. (To avoid confusion we note that Fig. 2 
is the distribution of total x 2 among 30 solutions, while 
Fig. 4 is the distribution of single-particle x 2 f° r the 29 
objects in the best solution for each case.) To reduce the 
spread of values we plot the square root of the sum of the 
squares of the normalized differences of model and catalog 
redshifts and distances. The curve is the distribution for 
Gaussian normal scatter of the redshifts and distances. 
There tend to be more objects with \ 2 < 1 than expected 
for this idealized model, even in the random case. This 
is a result of the relaxation to minimize x 2 , as discussed 
in connection with Figure 3. Even in the best solution 
for the real data there are more large values of x 2 than 
expected from the idealized distribution, but the excess is 
not large. For example, in cO there are 10 of the 29 objects 
with x 2 > 4, about 2.5 times the number expected from 
the idealized model. 

Our computation allows many attempts to find an ac- 
ceptable fit to catalog redshifts and distances, and the fit 
may be achieved at the expense of an unrealistic initial 
condition, associated with an unacceptable value or gra- 
dient of the primeval mass density fluctuation. Estimates 
of the primeval density fluctuations associated with a nu- 
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Fig. 4. — Distributions of single-particle \ 2 m redshift and dis- 
tance modulus for the 29 objects in the solutions with the smallest 
total x 2 f° r each case. In the broken histograms a = 0, in the solid, 
a = 15 km s — . 




Fig. 5. — Distributions of estimates of the local primeval density 
contijast extrapolated to the present in linear perturbation theory 
(eq. [|2C|]) in the best solutions for the six cases. The thin histograms 
assume u = 0, the thick, cu = 15 km s _1 . 



merical action solution are discussed in Peebles] (1996). 
We use a simplified approach based on the linear pertur- 
bation relation between the mass density contrast S and 
the peculiar velocity field v(x), 



dS 



-V- v. 



(19) 



This expression averaged over a sphere of comoving radius 
x is d = —3v r /ax, where v r is the peculiar velocity normal 
to and averaged over the surface of the sphere. The density 
contrast at high redshift varies as S cx D(t), where D(t) is 
the growing solution to the linear perturbation equation, 
so an estimate of the local density contrast extrapolated 
to the present in linear perturbation theory is 



So 



, Do (xj - ±j) ■ (xj - Xj) 

b{t) 4 



(20) 



The dot means derivative with respect to proper world 
time t, D a is the present value of D, and Xj and X; are 
the coordinate position and velocity of particle i at time 
t. We use the approximation x = (X3 — X2)/(*3 — £2), from 
the second to the third time step, because the fractional 
increase in time from the first to the second step is large, 
though it yields similar results. 

Figure 5 shows distributions of 5 from all pairs of ob- 
jects in the inner catalog for the best solution for each of 
the six cases. When the outer zone is removed from the 
dynamics the central value of the contrast 5 is close to 
unity, the value expected from the choice of M/L. The 
outer zone objects increase the central value of S for the 
inner zone, we suspect because the inner and outer zone 
objects are somewhat intermingled at high redshift. There 
are substantial differences among the distributions from 
solutions with similar values of x 2 ; reproducible features 



include the smaller central value of 5 when the outer zone 
is removed, the prominent positive tail in the distribution 
of 5 from the real catalogs, and the still more prominent 
tails in the random case. 

Because we have used a crude estimate of S a , and 
lumped together a range of values of primeval separations 
Xij, we could not expect the distributions in Figure 5 to 
be Gaussian even if we had an adequate approximation to 
Gaussian initial conditions. However, we might expect the 
more realistic cases to be closer to Gaussian. By this cri- 
terion the random control case is the least realistic, as we 
would hope. On the other hand, the outer zone increases 
the positive tail for the objects in the inner zone. This il- 
lustrates the difficulty of testing models at the present level 
of accuracy of the redshift and distance measurements. We 
turn now to a more demanding test. 

5. PROPER MOTIONS 



The satellite missions SIM flhttp:// sim.jpl.nasa.govj ) and 
G AIA ( http:/ /astro. estec.esa.nl/GAIA| ) could measure the 
velocities of nearby galaxies relative to the Milky Way and 
perpendicular to the line of sight, and thereby test the as- 
sumption that the galaxies trace the local mass distribu- 
tion. 

Figures 6 to 8 show examples of transverse velocities in 
our 30 solutions The larger symbols represent the better 
fits to the catalog redshifts and distances, with reduced \ 2 
for the inner objects less then 2.5, the smaller symbols the 
solutions with larger % 2 . The symbol orientations distin- 
guish solutions with and without the outer zone objects in 
the dynamics, as explained in the caption to Figure 6. 

When the dynamics includes only the inner catalog the 
transverse galactocentric velocity of M31 (Fig. 6) is close 
to zero in a few solutions, but in most cases M31 is moving 
at close to ±150 km s _1 normal to the supergalactic plane. 
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Fig. 6. — Transverse galactoccntric velocity of M31 in all 30 
solutions. The larger symbols have reduced total \ 2 l ess than 2.5. 
Squares, crosses and triangles pointing up use the two zones, case 
cl5. Diamonds, plus signs, and triangles pointing down show the 
effect of removing the outer zone. The symbol shapes are chosen 
according to «sgb, to facilitate comparison with the next 2 figures. 



This is part of a collapse to the local sheet-like arrange- 
ment shown in Figure 1. When the outer zone objects are 
included in the dynamics the fit is improved, the solutions 
have a broader range of values of the transverse velocity, 
and there are more solutions with a near radial approach 
of M31. The better solutions still prefer motion toward 
the supergalactic plane, and there are many possible di- 
rections of transverse motion that could appear to falsify 
the action model. 

The distribution of solutions for the transverse galac- 
tocentric velocity of IC 1613 in Figure 7 is broader when 
the outer zone is included, as for M31. One might have 
expected to see a correlation of solutions for the velocities 
of M31 and this distant companion of M31. There is no 
pronounced correlation in this example, but we regard this 
as a very preliminary indication to be considered in more 
detail. 

The orbit of NGC 6822 is more complex than M31 be- 
cause NGC 6822 is closer and moving away. Figure 8 shows 
one can find acceptable fits to the catalog redshift and 
distance, with a broader range of possible proper motions 
than in the previous two examples. We suspect this illus- 
trates a loss of memory of initial conditions in nonlinear 
classical dynamics. Tighter constraints on acceptable so- 
lutions might arise from a study of the initial conditions 
in the solutions; this also requires more study. 

6. DISCUSSION 

6.1. The Assumption Galaxies Trace Mass 

If the nearby galaxies trace the mass on scales ^> 100 kpc 
then our inner catalog in Table 1 gives a reasonably com- 
plete description of the mass within Mpc, and Ta- 
ble 2 gives a first approximation to the effect of external 



Fig. 7. — Transverse galactoccntric velocity of IC 1613. The 
symbol shapes are chosen according to the value of vgGB for M31, 
and the sizes according to the value of the total x 2 , as in Fig. 6. 
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Fig. 8. — Transverse galactocentric velocity of NGC 6822. The 
symbol shapes are chosen according to the value of vsgb f° r M31 
in Fig. 6, but here the larger symbols represent solutions in which 
the radial velocity of NGC 6822 is within 30 km s _1 of the catalog 
value. 



masses on the Local Group and its immediate neighbors. 
We have found dynamical solutions for the orbits of the as- 
sumed mass tracers that end up in our inner zone (within 
4/iy 5 Mpc distance); in the best cases the rms difference 
between model and catalog redshifts and distances is in 
the range 1.25 to 1.5 times the nominal catalog errors. 
The better fits use a more liberal estimate for the red- 
shift errors, an issue that requires further consideration. 
The difference of the reduced x 2 from unity is consider- 
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ably larger than statistical, but small enough to argue for 
the assumption galaxies trace mass. 

One of our controls keeps the catalog distances, red- 
shifts, and masses but places the objects at random posi- 
tions in the sky. Since this preserves the near Hubble flow 
of the more distant members of the inner zone and the 
infall of most of the closer members, it is not surprising 
that we get a rough fit. But this scrambling of the present 
angular positions worsens the fit, as expected if galaxies 
trace mass. This result is quite stable under changes in 
the random numbers for the scrambling and computation. 

Our second control removes the outer masses from the 
dynamics. It shows that the outer masses improve the fit 
to the measured redshifts and distances in the inner zone, 
again as expected if galaxies trace mass. This also suggests 
our solutions for the inner zone would be improved by 
using a better model for the external mass distribution, 
another point that requires further consideration. 

Figures 6 to 8 show that, if galaxies trace mass, there 
are solutions with similar values of \ 2 an d transverse ve- 
locities that Hiffer by 300 km s -1 These solutions tend 



to ocJnpy restricted regions of transverse velocity space 
The predictions for the galactocentric transverse velocity 
of M31 are only weakly perturbed by the masses beyond 
4:hj^ Mpc. Thus we suspect that under the assumption 
galaxies trace mass the allowed velocity of M31 is rather 
tightly constrained by the present catalogs. A measure- 
ment of the transverse velocity of this and other nearby 
galaxies could critically test the picture of the local mass 
distribution. 

6.2. Comments on Future Work 

This new version of the numerical action method could 
and c ertainly ehould bo applied to numerical N body aim 



ulatiobs — We hope to continue work on this topic — Our 



exploration of the effect of adjusting the mass-to-light ra- 
tios in the data has been quite limited, and here too we 
hope to be able to report on the values of M/L for late 
and early type systems that yield the best fit to catalog 
redshifts and distances, to compare to global measures of 
the density parameter. 

The allowance for errors in both redshift and distance 
is more realistic than previous applications of the action 
method, but it still ignores the errors in masses. A study of 
this issue certainly will be part of a satisfactory assessment 
of the dynamics of the nearby galaxies. 

A fundamental limitation of our numerical application 
should be noted. The matrix inverse of the second deriva- 
tives of the action with respect to the coordinates of an 
orbit gives the rate of change of the present distance with 
respect to the redshift (eq. ||]), and it reaches saddle points 
and local maxima as well as the local minima that are 
reached by following the gradient of the action. Because 
the computation time for the matrix inverse varies as the 
cube of the matrix size we relax the coordinates for each 
particle separately, toward a stationary point of the action 
and a minimum of the single-particle \ 2 ■ This is iterated 
through all particles many times. We are really interested 
in the minimum of the sum of the \ 2 fo r all particles, which 
need not minimize any single-particle x 2 - The inversion of 
the much larger matrix of all coordinates of all particles 
would be a heavy computation, but worth attempting. 

The galaxy NGC 6822 illustrates the problem of deal- 



ing with complex orbits. We have found reasonably close 
fits to its positive radial velocity by relaxing from ran- 
dom orbits, but it would be feasible and useful to develop 
a more systematic sample of the many allowed orbits for 
such cases. The results would challenge us to decide which 
solutions could have come from realistic initial conditions; 
we need something better than the simple estimator in 
equation (pp|). 

This work was supported in part at Princeton University 
by the NSF. 
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Table 1. Inner Catalog 
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Table 1. (Continued) Inner Catalog 


Name 


Note 


SGL 


SGB 


distance 


a v 


rcdshift 




mass 




Xz 
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Notes to Table 1. The following galaxies are combined into single entries. 

a) IC 342, NGC 1560, UGCA 105 

b) Maffei 1, Maffei 2 

c) M 31, M 33, IC 10, LGS 3 

d) M 81, M 82, NGC 2976, IC 2574 

e) NGC 5128 = Cen A, NGC 4945 

f) NGC 253, NGC 247 

g) NGC 5236 = M 83, NGC 5102 

h) NGC 2403, NGC 2366, Holmberg II 

i) NGC 55, ESO 294-010 
j) NGC 3109, Antlia dwarf 
k) Sextans A, Sextans B 
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Table 2. Outer Catalog 



Name Group SGL SGB distance ct m redshift a cz mass mass/dist a 
(deg) (deg) (Mpc) (km/s) (km/s) (log 10 [M Q ]) (1O 1O M Mpc' 3 ) 



Virgo 


11 


-1 


102. 


.7 


-2. 


.4 


17. 


.0 





.17 


986 


63 


15 


.1 


24 


Fornax 


51 


-1 


262 


.2 


-40. 


.9 


18. 


.8 





.20 


1466 


78 


14 


.5 


4 


Dorado 


53 


-1 


234 


.6 


-40 


.4 


18. 


.5 





.20 


1076 


65 


14 


.3 


3 


Coma I 


14 


-1 


87 


.0 


1. 


.2 


16. 


.4 





.40 


983 


53 


14 


.1 


3 


Ursa Major 


12 


-1 


65. 


.8 


2. 


.6 


16. 


.6 





.20 


957 


20 


13 


.9 


2 


Leo 


15 


-1 


93 


,7 


— 25. 


.8 


11 


.1 





.17 


72z 


61 


13 


.8 


4 


Sombrero 


11- 


14 


126. 


.2 


-6. 


.7 


10. 


.0 





.30 


985 


47 


13 


.5 


3 


N6744 


19 


-1 


207. 


.8 


10 


.2 


13 


.9 





.40 


752 


31 


13 


.3 


1 


CVn II 


14 


-4 


70. 


.5 


5 


.6 


7. 


7 





.20 


556 


12 


13 


.0 


2 


M51 


14 


-5 


73. 


,0 


16 


.1 


8. 


.0 





.40 


562 


46 


13 


.0 


2 


N1023 


17 


-1 


340 


.4 


-8. 


.0 


10. 


.0 





.20 


696 


10 


13 


.0 


1 


M101 


14 


-9 


63 


.3 


22. 


.6 


7. 


.4 





.20 


354 


27 


12 


.9 


2 


N6946 


14 


-0 


10. 


.0 


42. 


.0 


6. 








.60 


271 


6 


12 


.8 


3 


CVn I 


14 


-7 


76. 


.4 


6. 


.1 


4. 


9 





.30 


316 


11 


12 


.6 


4 



